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Doppler cooling the center-of-mass motion of an optically levitated microsphere via the velocity 
dependent scattering force from narrow whispering gallery mode (WGM) resonances is described. 
Light that is red detuned from the WGM resonance can be used to damp the center-of-mass motion 
in a process analogous to the Doppler cooling of atoms. Leakage of photons out of the microsphere 
when the incident field is near resonant with the narrow WGM resonance acts to damp the motion 
of the sphere. The scattering force is not limited by saturation, but can be controlled by the 
incident power. Cooling times on the order of seconds are calculated for a 20 micron diameter silica 
microsphere trapped within optical tweezers, with a Doppler temperature limit in the microKelvin 
regime. 

PACS numbers: 37.10.Vz, 37.30.+i, 42.50.Wk 



Doppler cooling has been an extremely successful tech- 
nique for cooling atomic species to temperatures in the 
microKelvin regime, opening up new areas in atomic 
molecular [2], condensed matter [sj, and many body 
physics It has allowed the creation of atomic gases 
in the quantum regime, including the creation of Bose- 
Einstein condensates of atomic gases [5] and Fermi gases 
[sl . More recently there has been considerable interest in 
the cavity cooling of atoms and molecules because a wider 
range of particles can, in principle, be cooled as no inter- 
nal resonance is required [7]. A resonance is, however, re- 
quired in the form of an external optical cavity, and a sin- 
gle atom [8], ion [9] and atomic ensembles [10] have been 
cooled. For molecular and atomic species that cannot be 
laser cooled, cavity cooling of a trapped species appears 
attractive because it does not rely on the detailed internal 
level structure. Over the last ten years the field of cavity 
optomechanics has cooled micro- and nano-mechanical 
objects down to temperatures where the quantum me- 
chanical nature of their motion will soon be apparent 
[lll-[l3] . Like cavity cooling of atoms and molecules, blue 
shifted photons are scattered from the cavity with respect 
to the incident photons, thus extracting energy. In this 
process, at least one degree of the freedom, such as a cav- 
ity mirror, or intracavity membrane, is d amp ed or cooled 
by interaction with the cavity field jisl, 18-21|. An im- 
portant system of this type is the cooling of the internal 
mechanical modes of a high Q, WGM resonator formed 
by a toroidal or spherical structure [22]. Very recently 
there have been proposals to cool optically levitated par- 
ticles using cavity cooling [23-25]. This is attractive be- 
cause the ability to levitate the particle isolates it well 
from the environment, increasing the prospects of cooling 
the center-of-mass motion to its quantum ground state. 
While this latter scheme is attractive for cooling nanopar- 
ticles, it does not appear to be practical for larger par- 
ticles, which would significantly perturb the cavity field 
reducing and potentially inhibiting cooling. 

In this letter we describe a hybrid scheme that links 
laser Doppler cooling and cavity cooling. This scheme 



differs from cavity cooling where the particle is cooled 
within the cavity, or in optomechanics, where part of the 
cavity is cooled. Instead, we cool the whole cavity in a 
process analogous to Doppler cooling where the required 
frequency dependent scattering force is provided by the 
high Q, WGM of the microsphere. 




FIG. 1: (Color online) The propagation of counter propagat- 
ing light rays (black and grey) from an incident plane wave 
within a microsphere on a whispering gallery mode resonance. 
The white rays are those from the incident field which im- 
pinge on the sphere but do not couple into the WGM. The 
corresponding isotropic leakage of photons from the sphere 
is analogous to the average isotropic spontaneous emission of 
photons from an atom following excitation. 

Whispering gallery modes or morphology dependent 
resonances occur in cylindrical and spherical dielectric 
particles which act as high Q (>10^) optical cavities for 
light that is coupled into the sphere and that propa- 
gates by total internal reflection around its annulus. Fig- 
ure 1 illustrates, from a geometrical optics perspective, 
the propagation of trapped rays of light from a incident 
plane wave within one plane through the microsphere. 
In spherical particles, the excitation of these WGM res- 
onances increases the scattering cross-section and there- 
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fore the radiation pressure that can be apphed to them. 
The highest Q's can be achieved when the radius of the 
microsphere far exceeds the wavelength of the hght cou- 
pled into it. The spectral widths of these resonances can 
vary considerably depending on wavelength and the size 
of the sphere but, like atomic resonances, they can have 
spectral widths of a few MHz, which are typically lim- 
ited by weak absorption or by Rayleigh scattering within 
the sphere. For these modes the sphere can be seen as 
a high Q spherical ring-like cavity. The effect of these 
resonances on the radiation pressure forces was first ob- 
served in the early work on levitating spheres by Ashkin 
and Dziedzic [26|. Although the coupling of light into 
these modes is not efficient for a free space propagating 
optical field, it was observed that the scattering cross- 
section, and also the optical force from levitation exper- 
iments, is enhanced when the incident light is resonant 
with these modes. Although not discussed by them in 
this context, the very narrow resonances could be used 
to damp the motion of the spheres via the Doppler effect, 
which transforms a frequency dependent force to a veloc- 
ity dependent force. It is stressed that it is the damping 
of the center-of-mass motion of the whole microsphere 
resonator structure and not the internal degrees of free- 
dom in spherical and toroidal resonators which has pre- 
viously been considered. 

To illustrate the forces on a resonator we first consider 
light incident on an idealized Fabry-Perot resonator, as 
shown in figure 2a, which is moving towards the beam. 
We consider the simple case where there are no losses in 
the mirrors or the medium between them and the reflec- 
tivity of each mirror is close to unity. The force on the 
cavity in a vacuum in the direction of propagation of an 
incident plane wave is calculated from the incident (Pi), 
reflected (Pr) and transmitted power (Pt), and is given 
by P = l/c(Pi — Pr — Pt). When the incident light is not 
near the cavity resonance almost all the incident light is 
reflected from the cavity and the force is at its maximum 
F ^ 2Pi/c. On resonance all light is transmitted and 
P ^ 0. The force is frequency dependent near the cavity 
resonance and is shown in figure 2a. The resonance is in- 
verted when compared to an atomic resonance, such that 
through the Doppler effect more force would be felt by 
the resonator when it is moving towards an incident field 
that is blue detuned with respect to the resonance. The 
WGM's in a microsphere are, however, more like the ring 
cavity as shown in single plane in figure 1 . Two rays near 
the annulus of the sphere are coupled and counter prop- 
agate around the sphere. To understand how this type 
of resonator is affected by radiation pressure we consider 
a more simplistic ring cavity of eight mirrors shown in 
figure 2b. Two rays of equal power from an incident field 
are coupled into the side mirrors. Any other rays from 
an optical field that are not resonant because of their 
angle of incidence or position on the cavity, will be re- 
flected and/or refracted producing a force that will not 



be strongly frequency dependent. We again consider the 
incident, reflected and transmitted light at each mirror, 
and the resulting forces due to the change in momentum, 
assuming that all mirrors have the same reflectivity and 
transmissivity with no losses. Unlike the Fabry-Perot 
resonator the transmitted rays are distributed evenly in 
a plane. By symmetry, the resulting forces around the 
ring due to the transmitted rays cancel each other out on 
average, just as in the case of isotropic spontaneous emis- 
sion from an excited atom. The force in the ^/-direction is 
then only due to the reflected and the incident field. Off 
resonance all the incident light is reflected and there is no 
net force in the ^/-direction. The force only acts to main- 
tain the sphere position along the x-axs. On resonance 
all light is transmitted and the net forces act to push 
the sphere in y direction, while maintaining the sphere 
on in its position on the x-axis. Figure 2b also shows 
a plot of the force in the ^/-direction due to these two 
rays derived from the reflection of an n mirrored cavity 
as function of frequency, where P = 2/c{P\ — Pr) • y and 

Pr = ^\(t — ^[ZlJ^-ikL )P and r is the amplitude re- 
flection coefficient [27]. In contrast to the Fabry-Perot 
cavity the force due to radiation pressure is maximised 
on resonance, just as in the atomic case, and the sphere 
acts in this respect like a large two level atom. 

The discussion above served to illustrate the basic 
physics behind the radiation pressure at resonance. We 
now calculate a more accurate force due to radiation pres- 
sure from an incident plane wave and the scattered light 
fields via Lorenz-Mie theory. For plane wave illumina- 
tion the radiation pressure on a non-absorbing sphere in 
vacuum, with size parameter x = ^a, where a is the 
sphere radius and A is the wavelength of light, is given 
by P = ^Qrad where P is the beam power and c is speed 
of light. The normalised radiation pressure cross-section 
is given by 

Qrad = Qext - , ' Re{anal^^ + 

^ n + 1 

n=0 

oo 

where Qext = 2/x^ ^(2n + l)Re{a 

n + ^n)- The val- 

n=0 

ues for the Mie coefficients and can be found from 
standard texts on Mie scattering [28j , where n represents 
the nth partial wave for the and hn modes respectively. 
The Mie resonances, or WGM's, for each can be found by 
solving for Im(an) = and Im(6n) = 0. The mode order / 
are the roots of the partial wave of mode number n, with 
the lowest order / = 1 producing the narrowest resonance 
for the each mode number. Figure 3 is a plot of the ra- 
diation pressure force calculated for an incident power of 
10 mW and size parameters from 39 to 41. A similar plot 
is obtained for a nearly collimated Gaussian beam using 
generalised Lorenz-Mie theory [29]. A range of resonances 
in the scattering force can be observed with very differ- 
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FIG. 2: (Color online) a) A diagram illustrating the incident 
(Pi), reflected (Pr) and transmitted (Pt) power from a Fabry- 
Perot cavity in a vacuum. Also shown is a diagram illustrating 
the variation in force due to radiation pressure as a function of 
frequency near the transmission peak. The vertical line rep- 
resents the detuning with respect to resonance, which would 
preferentially damp motion of the cavity towards the incident 
field, b) A diagram of an eight mirrored ring cavity to illus- 
trate the effect of radiation pressure force on WGM-like ring 
cavity which approximates the whispering gallery modes of a 
microsphere. Also shown is a diagram illustrating the varia- 
tion in force in the ^-direction as a function of frequency near 
resonance with the vertical line representing the red detuning 
required to damp the motion of the resonator moving in the 
-ve ^/-direction. The Lorentzian force profile, which peaks on 
resonance, is analogous to an atomic resonance in which the 
scattering force is peaked on resonance. 

ent widths. The narrowest are not well resolved in this 
low resolution calculation and are indicated by the solid 
vertical lines. Note that the size parameter can represent 
a variation in laser wavelength/frequency or a variation 
in particle radius a. The inset graph is a calculation at 
higher resolution on an expanded scale for the an( n = 52 
and / = 1) resonance for size parameter x = 40.62425. 
This would correspond to a 10 micron radius sphere illu- 
minated with light at approximately 773 nm. The scale 
has been converted to frequency with respect to the line 
centre of the resonance, which has a Lorentzian profile 
with a half width S = 27rx32 MHz. On resonance the 
force is 18.5 pN, which on this narrow frequency scale 
is offset by a constant force of 14.3 pN due to loight 
that is not reosnant with the WGM's. The force near 
any of these narrow Mie resonances can be approximated 

by Frad = ^ + ^ ioj-JoY'+s'' ' w^^^^ ^ constant 
over the frequency range considered and is the peak 
resonant force. Here 5 is the half- width- half-maximum 
line width and uj is the frequency of the light that is 
detuned from the resonant frequency ujq. The force is 



dependent on the velocity, of the microsphere via the 
Doppler effect and is given by Frad = ^ + ^-(A±£yTs^ 
where A is the detuning from resonance of the station- 
ary sphere. For the small velocities expected of a trapped 
microsphere, the force can be expanded about v = to 

give Prad = ^ + ^ + f (A^W""- ^^^^^^ ^^'^^ 

Doppler cooling of atoms there is no saturation of the 
cooling force. Therefore, in principle, the damping rate 
can be determined by the power of the light field Pp, 
although in practice additional heating, introduced by 
absorption, is likely to shift the resonance. 
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FIG. 3: (Color online) Calculated force on microsphere for 
size parameters from x = 39 - 41. This corresponds to a 
sphere with a radius of 10 /am illuminated by a plane wave 
of approximately 773 nm. Inset is a higher resolution of plot 
of the force as a function of frequency for the an{ n = 52 
and / = 1) resonance corresponding to a size parameter x = 
40.62425. 

We initially consider a single sphere cooled by a 1- 
D optical molasses using two counter-propagating fields. 
Here the force on the particle is P = where /3 = 
0^(1^+^^)^ . The e~^ velocity damping time or cooling 
time is r ^ {/3/m)~^. Using Po=100 mW, beam and 
mass m = 4 X 10~^^ kg corresponding to a =10 /im Si02 
sphere, gives a characteristic cooling time of r = 5.7 
s for (5 = 27rx32 MHz HWHM resonance or 178 ms 
for a 27rxl MHz resonance width. Like laser Doppler 
cooling of atoms, this will lead to a 1-D cooling limit 
based on a balance between the average energy damp- 
ing rate or cooling power expressed as (P) = f3{v'^) and 
heating by diffusion with the diffusion constant given by 
D = {p^) = h^k'^Vsc where Vsc is the frequency depen- 
dent scattering rate given by Fgc = ^ {z^+Ly^+s'^ - 
where l/2ki,T=l/2m{v'^) and where A = (5, the Doppler 
cooling limit is the same as for atomic systems and is 
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given by ^ hS. For the 32 MHz resonance this cor- 
responds to a temperature of 760 /iK, and 24 /iK for the 
1 MHz resonance. 

An efficient optical trap for dielectric microspheres is 
produced by two diverging counter-propagating fields in 
a two-beam fiber trap [30|. A single collimated beam 
that is red detuned from resonance could be used to cool 
the trapped sphere in analogy to the Doppler cooling of a 
trapped ion. As microspheres are often trapped in air we 
compare the damping of their motion due to gas viscosity 
and to Doppler cooling. For an optically sphere trapped 
within a viscous medium such as air, the equation of 
motion along one axis of the trap is given by 

mx'\t) = -ujlx{t) + (A - To)x\t) + Fo + Ff (t), (1) 

where cuq = ^J nlm is the trap frequency, k is the spring 
constant of the trap. The drag coefficient for damping of 
a spherical particle by gas viscosity {r\) is Fq = 67r77a and 
by Doppler cooling is /3t = ^(A^+^^p • The microsphere 
will also be subjected to a frequency independent offset 
force Fo = + a^+c^^ ^ time- varying Langevin 
force Ff{t). At 288 K, the damping in air is given by F = 
3.4 X 10~^ kgs~^ while the maximum value of the optical 
damping coefficient for a 100 mW incident beam is only 
three orders of magnitude less at /3t = 6.7 x 10""*^^ kgs~^. 
These values are approximately equal at a pressure of 15 
mTorr, based on the drag coefficient of a sphere in the free 
molecular fiow regime where Fq = (4/3 + 37r/16)7rp('u)a^, 
{v) is the mean velocity of the gas particles and p is the 
gas density 0. The optical damping is six orders of 
magnitude greater than that due to the background gas 
at 1x10-^ torr. 

A typical spring constant of an optical trap is k = 
5xl0~^ and uJo=2 ttx 740 Hz and, using the mass cal- 
culated from a spherical silica sphere of radius a = 10 
microns, the motion should be damped on the order of 
seconds. This does not correspond to sideband resolved 
cooling, but this could be accomplished by cooling on a 
resonance that is narrower than the trap frequency or 
by increasing the trap frequency via an increase in light 
intensity used to trap the particle. The equation of mo- 
tion above also applies to a microsphere attached to a 
cantilever and Doppler cooling of the microsphere can 
be used to cool the cantilever motion. Here the spring 
constant of the cantilever is much larger than an optical 
trap (<80 N/m) with oscillation frequencies below 1 MHz 
typically in the 100 kHz range. A cantilever of A: = 77 
N/m and resonant frequency of 1 MHz, has an effective 
mass of 2 X 10~^^ kg which is similar to the mass of the 
10 fim sphere considered above. Sideband resolved cool- 
ing therefore also appears feasible in this optomechanical 
system. 

As the resonance frequency for any microsphere is criti- 
cally dependent on its radius, the resonance condition for 
each sphere must be found in order to begin cooling it. 



This requires that the laser must be offset locked to the 
resonance. Locking a laser to WGM resonance has been 
demonstrated [32[ and this will be more easily found for 
the microsphere cantilever system since the sphere is per- 
manently attached to the cantilever. Silica microspheres 
have been trapped in vacuum for up to half an hour in 
vacuum at 10 ~^ torr with the time limited by radiomet- 
ric heating from absorption of the laser light at 514.5 nm 
(33} . This time could be considerably improved by us- 
ing low loss Si02 spheres, as well as using trapping and 
cooling light in the low-loss wavelength window around 
1-1.5 microns where narrow bandwidth low noise lasers 
are available. Finding and subsequently locking the laser 
to the appropriate WGM could be carried out while the 
particle is trapped in air, which we have demonstrated 
can be trapped in excess of 5 hours. Once locked to the 
red side of the resonance, the air could be pumped out 
and the microsphere cooled to its Doppler limit. A micro- 
sphere has internal mechanical resonances which can be 
excited by radiation pressure acting on the internal sur- 
faces of the sphere when light is coupled in the WGmQ. 
This will induce motional sidebands on the microspheres 
WGM resonance. When light is red detuned from res- 
onance it has been shown that this motion can also be 
cooled, and thus it may be feasible to cool both the in- 
ternal and external degrees of freedom the microsphere. 
This may help to ameliorate the very small residual ab- 
sorption of light which will act to heat the sphere. 

A method for cooling the center-of-mass motion of a 
microsphere using the velocity dependent force inherent 
in the whispering gallery modes was presented. This type 
of Doppler cooling has much in common with laser cool- 
ing of trapped atoms and ions. A 3-D optical molasses 
could be used to cool a microsphere in all three dimension 
or by a single beam when the sphere is trapped using opti- 
cal or electrostatic fields or attached to a cantilever. Such 
a scheme may also be used to sympathetically cool op- 
tically bound, co-trapped particles which do not possess 
whispering gallery mode resonances. Like laser cooling 
ultimate temperatures in the microkelvin range appear 
feasible with cooling times on the order of seconds. 
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